GYRO.013A PATENT 
SYSTEM AND METHOD FOR MEASUREMENTS OF DEPTH AND VELOCITY OF 
INSTRUMENTATION WITHIN A WELLBORE 



Claim of Priority 

[0001] This application claims benefit to U.S. Provisional Application No. 
60/539,234, filed January 26, 2004, which is incorporated in its entirety by reference herein. 

Background of the Invention 

Field of the Invention 

[0002] The present application relates generally to systems and method for 
determining the depth, the velocity, or both the depth and the velocity of an instrumentation 
package within a wellbore. 
Description of the Related Art 

[0003] Surface-based wellbore depth measurements are typically made 
periodically during wellbore drilling for the exploration of oil and gas deposits to determine 
the absolute depth of the drilling tool within the wellbore. In such measurements, the depth 
of the drilling tool is typically determined by surface measurements of the lengths of the pipe 
sections inserted into the wellbore between the drilling tool and the surface. 

[0004] Using wireline surveys, the drilling of the wellbore is periodically halted 
and a survey tool is lowered into the wellbore. As the survey tool is guided along the 
wellbore, it can provide information regarding its orientation and location by sending signals 
through a wire or cable to the surface. The absolute depth of the survey tool down the 
wellbore is typically given by a surface-based measurement of the length of the wire or cable 
between the survey tool and the surface. Similarly, surface-based logging measurements of 
the absolute depth of detected geological formations are typically made by surface 
measurements of the length of the wire or cable between the logging tool and the surface. 
However, due to various distortions, the cumulative length of the pipe sections or of the wire 
while within the wellbore can differ from the cumulative length measured at the surface, 
resulting in errors in the determination of the absolute depth. 

[0005] In addition, surface-based measurements of the velocity of the survey tool 
along the wellbore can be used to determine the relative distances between detected features 
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or formations within the wellbore. For example, the time period between detecting two 
separate features along the wellbore and the velocity of the survey tool during this time 
period can be multiplied together to provide the relative distance between the two detected 
features. However, as with the surface-based depth measurements described above, surface- 
based velocity measurements do not provide a sufficient accuracy (e.g., within only a few 
centimeters) to tell when two geophysical sensors with significant along-hole separation pass 
the same geological formation within only a few centimeters of accuracy. For example, due 
to friction or other effects within the wellbore, the survey tool can move in a jerking manner 
with varying velocity. In addition, these effects can result in the velocity of the survey tool 
within the wellbore differing from the measured velocity of the wire or cable at the surface. 

[0006] Inertial navigation systems have been proposed as being able to provide 
improved wellbore depth and velocity measurements. Such inertial navigation systems can 
determine the depth of the survey tool by double integration of the detected acceleration. 
However, such procedures are vulnerable to noise in the detected acceleration, which results 
in a drift of the depth measurement from the true depth of the survey tool. It has been 
difficult, and will likely continue to be difficult, to get such inertial navigation systems to 
work without providing aiding data such as surface-based depth measurements or updates of 
the depth obtained while the survey tool is stationary (i.e., zero-velocity updates) to remove 
this drift. However, it is generally desirable to avoid using surface-based aiding data, since, 
such data would effectively transform an inertial navigation system into a surface-depth 
system, thereby including the problems of such systems. It is also generally desirable to 
avoid zero-velocity updates due to large increases in the time and costs associated with such 
surveys. 

[0007] Improved wellbore depth and velocity measurements are desirable to 
understand the geological formations being drilled and the oil or gas deposits being accessed. 
For example, improved measurements can resolve uncertainties in the depth measurements of 
a geological fault from wellbore positioning surveys in two nearby wells. In addition, 
improved wellbore depth measurements are helpful for drilling safety by providing more 
reliable information regarding the true wellbore depth to avoid drilling into adjacent wells. 
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Summary of the Invention 

[0008] In certain embodiments, a survey tool is adapted for use in a wellbore. 
The survey tool comprises a downhole portion having an axis. The downhole portion is 
adapted to move within the wellbore with the axis generally parallel to the wellbore. The 
survey tool further comprises a first acceleration sensor mounted at a first position within the 
downhole portion. The first acceleration sensor is adapted to generate a first signal indicative 
of an acceleration of the first acceleration sensor along the axis. The survey tool further 
comprises a second acceleration sensor mounted at a second position within the downhole 
portion. The second position is spaced from the first position by a non-zero distance along 
the axis. The second acceleration sensor is adapted to generate a second signal indicative of 
an acceleration of the second acceleration sensor along the axis. The survey tool further 
comprises a controller adapted to receive the first signal and the second signal and to 
calculate a depth, a velocity, or both a depth and a velocity of the downhole portion in 
response to the first signal and the second signal. 

[0009] In certain embodiments, a method determines a depth of a downhole 
portion of a survey tool along a wellbore. The method comprises providing a survey tool 
comprising a downhole portion. The downhole portion comprises a first acceleration sensor 
and a second acceleration sensor. The first acceleration sensor is adapted to generate a first 
signal indicative of an acceleration of the first acceleration sensor along the wellbore. The, 
second acceleration sensor is adapted to generate a second signal indicative of an acceleration 
of the second acceleration sensor along the wellbore. The second acceleration sensor is 
spaced from the first acceleration sensor by a non-zero distance. The method further 
comprises receiving the first signal and the second signal while the downhole portion is at a 
first location within the wellbore. The method further comprises receiving the first signal 
and the second signal while the downhole portion is at a second location within the wellbore. 
The method further comprises calculating a depth of the downhole portion of the survey tool 
in response to the first signal and the second signal received while the downhole portion is at 
the first location and in response to the first signal and the second signal received while the 
downhole portion is at the second location. 



[0010] In certain embodiments, a method determines a velocity of a downhole 
portion of a survey tool along a wellbore. The method comprises providing a survey tool 
comprising a downhole portion. The downhole portion comprises a first acceleration sensor 
and a second acceleration sensor. The first acceleration sensor is adapted to generate a first 
signal indicative of an acceleration of the first acceleration sensor along the wellbore. The 
second acceleration sensor is adapted to generate a second signal indicative of an acceleration 
of the second acceleration sensor along the wellbore. The second acceleration sensor is 
spaced from the first acceleration sensor by a non-zero distance. The method further 
comprises receiving the first signal and the second signal while the downhole portion is at a 
first location within the wellbore. The method further comprises receiving the first signal 
and the second signal while the downhole portion is at a second location within the wellbore. 
The method further comprises calculating a velocity of the downhole portion of the. survey 
tool in response to the first signal and the second signal received while the downhole portion 
is at the first location and in response to the first signal and the second signal received while 
the downhole portion is at the second location. 

[0011] For purposes of summarizing the invention, certain aspects, advantages and 
novel features of the invention have been described herein above. It is to be understood, 
however, that not necessarily all such advantages may be achieved in accordance with any 
particular embodiment of the invention. Thus, the invention may be embodied or carried out in 
a manner that achieves or optimizes one advantage or group of advantages as taught herein 
without necessarily achieving other advantages as may be taught or suggested herein. 

Brief Description of the Drawings 

[0012] Figure 1 schematically illustrates a survey tool compatible with 
embodiments described herein for use in a wellbore. 

[0013] Figure 2 schematically illustrates the survey tool in a portion of the 
wellbore having a curvature. 

[0014] Figure 3 schematically illustrates a survey tool as part of a logging 
assembly having a first supplementary sensor and a second supplementary sensor. 

[0015] Figure 4A is a flowchart of an exemplary method of determining a depth 
of a downhole portion of a survey tool in accordance with embodiments described herein. 
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[0016] Figure 4B is a flowchart of an exemplary method of determining a velocity 
of a downhole portion of a survey tool in accordance with embodiments described herein. 

[0017] Figure 5 A schematically illustrates a survey tool having a first 
supplementary sensor and a second supplementary sensor passing a landmark feature. 

[0018] Figure 5B is a plot of the signals from the first supplementary sensor and 
the second supplementary sensor as a function of time. 

Detailed Description of the Preferred Embodiment 

[0019] Certain embodiments described herein provide a true downhole-based 
system for measuring a depth, a velocity, or both a depth and a velocity of a downhole 
portion with sufficient accuracy for logging and drilling applications. 

[0020] Figure 1 schematically illustrates a survey tool 10 compatible with 
embodiments described herein for use in a wellbore 20. The survey tool 10 comprises a 
downhole portion 30 having an axis 32. The downhole portion 30 is adapted to move within 
the wellbore 20 with the axis 32 generally parallel to the wellbore 20. The survey tool 10 
further comprises a first acceleration sensor 40 mounted at a first position 42 within the 
downhole portion 30. The first acceleration sensor 40 is adapted to generate a first signal 
indicative of an acceleration of the first acceleration sensor 40 along the axis 32. The survey 
tool 10 further comprises a second acceleration sensor 50 mounted at a second position 52 
within the downhole portion 30. The second position 52 is spaced from the first position 42 
by a non-zero distance B along the axis 32. The second acceleration sensor 50 is adapted to 
generate a second signal indicative of an acceleration of the second acceleration sensor 50 
along the axis 32. The survey tool 10 further comprises a controller 60 adapted to receive the 
first signal and the second signal. The controller 60 is further adapted to calculate a depth, a 
velocity, or both a depth and a velocity of the downhole portion 30 in response to the first 
signal and the second signal. In the embodiment schematically illustrated by Figure 1, the 
controller 60 is at the surface and is coupled to the downhole portion 30 by a cable 62. 

[0021] In certain embodiments, the survey tool 10 is a component of a drill string 
and is used to determine the actual depth of a drilling tool (e.g., drill bit) of the drilling 
assembly. Drill strings compatible with embodiments described herein include, but are not 
limited to, measurement-while-drilling (MWD) strings. In certain other embodiments, the 
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survey tool 10 is a component of a navigational string and is used to determine at least a 
portion of the wellbore path. In certain other embodiments, the survey tool 10 is a 
component of a logging string and is used to determine the actual depth of detected 
geological features along the wellbore 20 or relative depths between detected geological 
features along the wellbore 20. Logging strings compatible with embodiments described 
herein include, but are not limited to logging-while-drilling (LWD) strings. In certain 
embodiments, the drill string or the logging string includes a sufficient number of sensors and 
adequate spacings between the first acceleration sensor 40 and the second acceleration sensor 
50 to perform the method described below. In certain embodiments, the drill string or the 
logging string includes redundant acceleration sensors (e.g., cross-axial accelerometers). 

[0022] In certain embodiments, the downhole portion 30 comprises a housing 
containing at least one of the acceleration sensors. As schematically illustrated by Figure 1, 
the housing of certain embodiments contains both the first acceleration sensor 40 and the 
second acceleration sensor 50. In other embodiments, the first acceleration sensor 40 and the 
second acceleration sensor 50 are not contained in a single housing, but are positioned on 
different portions of the downhole portion 30. 

[0023] In certain embodiments, the downhole portion 30 is adapted to bend as the 
downhole portion 30 moves through a curved portion of the wellbore 20. Figure 2 
schematically illustrates a downhole portion 30 within a section of the wellbore 20 having a 
curvature such that the direction of the wellbore 20 changes by a non-zero angle 0. In such 
embodiments, the downhole portion 30 bends by the non-zero angle #such that the axis 32 of 
the downhole portion 30 is substantially parallel to the wellbore 20. Under such conditions, 
the axis 32 at the first position 42 is at the non-zero angle with respect to the axis 32 at the 
second position 52. 

[0024] In certain embodiments, the first acceleration sensor 40 and the second 
acceleration sensor 50 comprise accelerometers currently used in conventional wellbore 
survey tools. In certain embodiments, one or both of the first acceleration sensor 40 and the 
second acceleration sensor 50 comprise a single-axis accelerometer sensitive to accelerations 
along a single sensing direction. In such embodiments, the single-axis accelerometer is 
advantageously mounted so that its sensing direction is substantially parallel with the axis 32 
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of the downhole portion 30. In other embodiments, one or both of the first acceleration 
sensor 40 and the second acceleration sensor 50 comprise a two-axis (i.e., cross-axial) 
accelerometer sensitive to accelerations in an x-y plane. In such embodiments, the two-axis 
accelerometer is advantageously mounted so that its x-y plane is parallel to the axis 32 of the 
downhole portion 30. In certain embodiments, the first acceleration sensor 40 and the second 
acceleration sensor 50 are advantageously substantially identical. Exemplary accelerometers 
include, but are not limited to, quartz flexure suspension accelerometers available from a 
variety of vendors. Other types of acceleration sensors are also compatible with 
embodiments described herein. 

[0025] In certain embodiments, the distance B between the first position 42 and 
the second position 52 along the axis 32 of the downhole portion 30 is advantageously 
selected to be long enough to respond to curvature of the wellbore 20 such that the downhole 
portion 30 bends substantially equally to the curvature of the wellbore 20. The distance B of 
certain embodiments is selected such that the first acceleration sensor 40 and the second 
acceleration sensor 50 detect different accelerations along the axis 32 when the downhole 
portion 30 is in a curved portion of the wellbore 20. In certain embodiments, the distance B 
is larger than approximately 10 meters, while in other embodiments, the distance B is in a 
range between approximately 10 meters and approximately 30 meters. Other distances B are 
compatible with embodiments described herein. 

[0026] In certain embodiments, the downhole portion 30 comprises one or more 
supplementary sensors in addition to the first acceleration sensor 40 and the second 
acceleration sensor 50. Figure 3 schematically illustrates a survey tool 10 comprising a first 
supplementary sensor 70 and a second supplementary sensor 80. Exemplary supplementary 
sensors of the downhole portion 30 include, but are not limited to, gamma-ray sensors 
adapted to detect gamma rays from geological formations in proximity to the downhole 
portion 30 and magnetic sensors adapted to detect casing collars of the pipe casing sections of 
the wellbore 20 in proximity to the downhole portion 30. In the embodiment schematically 
illustrated by Figure 3, the first supplementary sensor 70 is below and in proximity to the first 
acceleration sensor 40 and the second supplementary sensor 80 is above and in proximity to 
the second acceleration sensor 50. Other embodiments can have other configurations of at 
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least one supplementary sensor and the two acceleration sensors. In certain embodiments, 
such supplementary sensors are used in conjunction with the survey tool 10, as described 
more fully below, to provide additional data which is used to aid the determination of the 
depth, the velocity, or both the depth and the velocity of the downhole portion 30 within the 
wellbore 20. 

[0027] In certain embodiments, the controller 60 is adapted to determine the 
depth, the velocity, or both the depth and the velocity of the downhole portion 30 in response 
to signals received from the various sensors of the downhole portion 30. In certain 
embodiments, the controller 60 comprises a microprocessor adapted to perform the method 
described below for determining the depth, the velocity, or both the depth and the velocity of 
the downhole portion 30 within the wellbore 20. In certain embodiments, the controller 60 
further comprises a memory subsystem adapted to store at least a portion of the data obtained 
from the various sensors. The controller 60 can comprise hardware, software, or a 
combination of both hardware and software to accomplish the calculation of the depth, the 
velocity, or both the depth and the velocity. In certain embodiments, the controller 60 
comprises a standard personal computer. 

[0028] In certain embodiments, at least a portion of the controller 60 is located 
within the downhole portion 30. In certain other embodiments, at least a portion of the 
controller 60 is located at the surface and is coupled to the downhole portion 30 within the 
wellbore 20 by a wire or cable 62. In such embodiments, the cable 62 comprises signal 
conduits through which signals are transmitted from the various sensors within the downhole 
portion 30 to the controller 60. In certain embodiments in which the controller 60 is adapted 
to generate control signals for the various components of the survey tool 10 in the downhole 
portion 30, the cable 62 is adapted to transmit the control signals from the controller 60 to the 
downhole portion 30. 

[0029] In certain embodiments, the controller 60 is adapted to perform a post- 
processing analysis of the data obtained from the various sensors of the survey tool 10. In 
such post-processing embodiments, data is obtained and saved from the various sensors of 
the survey tool 10 as the downhole portion 30 travels within the wellbore 20, and the saved 
data are later analyzed to determine relative depths and/or absolute depths of the various 



detected features. The saved data obtained from the various sensors, including any aiding 
data (described more fully below) advantageously includes time reference information (e.g., 
time tagging) so that the relative times of the detection of various features can be determined. 
In certain embodiments, the relevant data from the various sensors are manually inspected 
and correlated with one another to provide aiding data. In other embodiments, the controller 
60 performs this correlation of the saved data automatically to provide the aiding data. 

[0030] In certain other embodiments, the controller 60 provides a real-time 
processing analysis of the data obtained from the various sensors of the survey tool 10. In 
such real-time processing embodiments, data obtained from the various sensors of the survey 
tool 10 are analyzed while the downhole portion 30 travels within the wellbore 20. In certain 
embodiments, at least a portion of the data obtained from the various sensors is saved in 
memory for analysis by the controller 60. The controller 60 of certain such embodiments 
comprises sufficient data processing and data storage capacity to perform the real-time 
analysis. In certain embodiments, the relevant data from the various sensors are 
advantageously correlated with one another by the controller 60 to provide aiding data. 

[0031] Figure 4 A is a flowchart of an exemplary method 100 for determining a 
depth of a downhole portion 30 of a survey tool 10 along a wellbore 20 in accordance with 
embodiments described herein. While the method 100 is described herein by reference to the 
survey tool 10 schematically illustrated by Figures 1-3, other survey tools 10 are also 
compatible with embodiments of the method 100. 

[0032] In certain embodiments, the method 100 comprises providing the survey 
tool 10 comprising a downhole portion 30 within the wellbore 20 in an operational block 
110. The downhole portion 30 comprises a first acceleration sensor 40 and a second 
acceleration sensor 50. The first acceleration sensor 40 is adapted to generate a first signal 
indicative of an acceleration of the first acceleration sensor 40 along the wellbore 20. The 
second acceleration sensor 50 is adapted to generate a second signal indicative of an 
acceleration of the second acceleration sensor 50 along the wellbore 20. The second 
acceleration sensor 50 is spaced from the first acceleration sensor 40 by a non-zero distance. 

[0033] In certain embodiments, the method 100 further comprises receiving the 
first signal and the second signal while the downhole portion 30 is at a first location within 
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the wellbore 20 in an operational block 120. The downhole portion 30 of certain 
embodiments is stationary while at the first location, while in other embodiments, the 
downhole portion 30 is moving along the wellbore 20 while at the first location. 

[0034] In certain embodiments, the method 100 further comprises receiving the 
first signal and the second signal while the downhole portion 30 is at a second location within 
the wellbore 20 in an operational block 130. The downhole portion 30 of certain 
embodiments is stationary while at the second location, while in other embodiments, the 
downhole portion 30 is moving along the wellbore 20 while at the second location. 

[0035] In certain embodiments, the method 100 further comprises calculating a 
depth of the downhole portion 30 in an operational block 140. The depth is calculated in 
response to the first signal and the second signal received while the downhole portion 30 is at 
the first location and in response to the first signal and the second signal received while the 
downhole portion 30 is at the second location. 

[0036] Figure 4B is a flowchart of an exemplary method 200 for determining a 
velocity of a downhole portion 30 of a survey tool 10 between two locations along a wellbore 
20 in accordance with embodiments described herein. While the method 200 is described 
herein by reference to the survey tool 10 schematically illustrated by Figures 1-3, other 
survey tools 10 are also compatible with embodiments of the method 200. 

[0037] In certain embodiments, the method 200 comprises providing the survey 
tool 10 comprising a downhole portion 30 in an operational block 210. The downhole 
portion 30 comprises a first acceleration sensor 40 and a second acceleration sensor 50. The 
first acceleration sensor 40 is adapted to generate a first signal indicative of an acceleration of 
the first acceleration sensor 40 along the wellbore. The second acceleration sensor 50 is 
adapted to generate a second signal indicative of an acceleration of the second acceleration 
sensor 50 along the wellbore. The second acceleration sensor 50 is spaced from the first 
acceleration sensor 40 by a non-zero distance. 

[0038] In certain embodiments, the method 200 further comprises receiving the 
first signal and the second signal while the downhole portion 30 is at a first location within 
the wellbore 20 in an operational block 220. The downhole portion 30 of certain 
embodiments is moving along the wellbore 20 while at the first location. In certain 
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embodiments, the method 200 further comprises receiving the first signal and the second 
signal while the downhole portion 30 is at a second location within the wellbore 20 in an 
operational block 230. The downhole portion 30 of certain embodiments is moving along the 
wellbore 20 while at the second location. 

[0039] In certain embodiments, the method 200 further comprises calculating a 
velocity of the downhole portion 30 between the first location and the second location in an 
operational block 240. The velocity is calculated in response to the first signal and the 
second signal received while the downhole portion 30 is at the first location and in response 
to the first signal and the second signal received while the downhole portion 30 is at the 
second location. 

[0040] An exemplary embodiment for determining the depth, the velocity, or both 
the depth and the velocity of a downhole portion 30 of a survey tool 10 utilizing a first 
acceleration sensor 40 and a second acceleration sensor 50 is described below. While the 
exemplary embodiment described below has a minimum number of variables, other 
embodiments are not limited to only these variables. Additional variables may also be used, 
including, but not limited to, misalignments of the acceleration sensors relative to the axis 32. 
In certain embodiments, the units of the parameters and variables below are in meters- 
kilogram-second (MKS) units. 

[0041] Aiding data from other downhole supplementary sensors (e.g., gamma-ray 
sensors, magnetic sensors for locating casing collars) is advantageously included in certain 
embodiments to enhance the resultant accuracy of the results. Other embodiments do not 
utilize aiding data or such supplementary sensors. The exemplary embodiment described 
below includes the use of aiding data, such as velocity data, absolute-depth data, and relative- 
depth data. Other types and/or combinations of aiding data are also used in other 
embodiments. 

[0042] In the exemplary embodiment described below, the periodicity of the 
measurements from the two accelerometers define time periods or "epochs" whereby one set 
of accelerometer measurements are taken at every epoch k. Aiding data are taken in the 
exemplary embodiments only at a subset of these epochs. In certain embodiments, the 
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different types of aiding data are taken the same epochs, while in other embodiments, the 
different types of aiding data are taken at different epochs. 

[0043] In the exemplary embodiment described below, the first acceleration 
sensor 40 is referred to as the "upper acceleration sensor" and the second acceleration sensor 
50 is referred to as the "lower acceleration sensor." The terms "upper" and "lower" are used 
herein merely to distinguish the two acceleration sensors according to their relative positions 
along the wellbore 20, and are not to be interpreted as limiting. Other embodiments 
distinguish the two acceleration sensors from one another using other terms. 
Exemplary Embodiment: State Vector With Five Elements 

[0044] The exemplary embodiment described below utilizes a state vector having 

five elements and the following parameters: 

g = magnitude of gravity; 

At = time between updates; and 

B = distance between the upper acceleration sensor 40 (denoted below by the 
subscript "U") and the lower acceleration sensor 50 (denoted below by the 
subscript "L"). 

In certain embodiments, the time between updates At corresponds to the clock frequency of 
the computer system used to perform the exemplary embodiment. 

[0045] The state vector^ at epoch k is expressed as follows: 
- **=k D Lk d k I Lk ] T ; (Eq.l) 
where a* is the calculated acceleration of the survey tool 10 in a direction generally parallel to 
the wellbore 20, is the calculated velocity of the survey tool 10 in a direction generally 
parallel to the wellbore 20, D L ,* is the calculated depth of the lower acceleration sensor 50, dk 
is the calculated apparent dogleg, which is equal to the difference between the inclinations of 
the lower acceleration sensor 50 and the upper acceleration sensor 40, divided by the distance 
B (i.e., d k -(l Lk -I Uk )/B), and I Lt k is the calculated inclination of the lower acceleration 

sensor 50. 
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[0046] The state co- variance matrix at epoch k is expressed as follows: 
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(Eq. 2) 



where o^* is the variance of parameter number i in state vector A*, and a iJt k is the co-variance 
between parameter number i and j in state vector A*. 

[0047] The -initial state at epoch = 0, corresponding to the survey tool 10 in a 
stationary condition within the wellbore 20 is given by the following: 

* 0 =[0 0 D Lfi (l Lfi -I Ufi )/B I L J; (Eq.3) 

where D Lt o is the initial depth of the lower acceleration sensor 50, Il,o is the initial inclination 
of the lower acceleration sensor 50, and I Ut o is the initial inclination of the upper acceleration 
sensor 40. 

[0048] In certain embodiments, the initial depth is referred to a known point in the 
wellbore 20. The initial inclinations of certain embodiments is determined from stationary 
acceleration sensor measurements or from known wellbore geometry of a landmark location. 
Additional acceleration sensors are compatible with the use of acceleration-sensor-based 
initial inclinations. In certain embodiments, a pair of high-side cross-axial acceleration 
sensors can be used to provide initial stationary acceleration measurements in a direction 
substantially perpendicular to the axis 32 and having a component in a vertical plane. In 
other embodiments, the first acceleration sensor 40 and the second acceleration sensor 50 
each comprises a two-axis accelerometer that provides signals indicative of the acceleration 
parallel to the axis 32 and the acceleration in a direction substantially perpendicular to the 
axis 32 and having a component in a vertical plane. In certain such embodiments, the initial 
inclinations are then given by: 

I L 0 = arctan(- HA L Q /A L 0 ); (Eq. 4) 

I uo = arctan(- HA U0 /A U 0 ); (Eq. 5) 

where Al,o is the initial lower acceleration sensor measurement in a direction generally 
parallel to the axis 32, A Ut o is the initial upper acceleration sensor measurement in a direction 
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generally parallel to the axis 32, HA L 0 is the initial lower acceleration sensor measurement in 
a direction substantially perpendicular to the axis 32 and having a component in a vertical 
plane (i.e., the lower high-side acceleration), and HAu,o is the initial upper acceleration sensor 
measurement in a direction substantially perpendicular to the axis 32 and having a component 
in a vertical plane (i.e., the upper high-side acceleration). 

[0049] The co-variance matrix Z 0 for the initial state at epoch k = 0 can be 
expressed as the following: 
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(Eq. 6) 



where ao is the uncertainty in the initial depth of the lower acceleration sensor 50 and <j\ is 
the uncertainty in the initial inclination of the lower acceleration sensor 50. The zero 
elements of the co-variance matrix Zq result from the fact that the downhole portion 30 is 
initially stationary (i.e., acceleration and velocity both equal zero). 

[0050] The state vector X k -i of epoch k-1 can be used to predict the state vector X k 
of a later epoch k using the following equations: 

(Eq. 7) 
(Eq. 8) 
(Eq. 9) 
(Eq. 10) 
(Eq.ll) 

where At is the time period between epoch k-1 and epoch k. In addition, other embodiments 
utilize other equations which include higher-order terms to predict the state vector of epoch k 
based on an earlier state vector of epoch k-1. 

[0051] The co-variance matrix % for the predicted state vector is given by the 
following diagonal matrix: 
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where p a is the maximum change of acceleration over time period At, p v is the maximum 
change of velocity over time period At, po is the maximum change of depth over time period 
At, pd is the maximum change of apparent dogleg over time period At, and pi is the maximum 
change of inclination over time period At. In certain embodiments, p„ is assumed to be given 
by p a = 2po/(At) 2 . In certain embodiments, p v is assumed to be given by p v = polAt. In 
certain embodiments, pa is assumed to be given by pd = 2p//B. 

[0052] The parameter a provides a multiplication factor between the standard 
deviation crof a state vector element and the maximum change p of the state vector element, 
such that the maximum change of the state vector element can be expressed as p = a* a. In 
certain embodiments, the multiplication factor or is in a range between approximately 2 and 
approximately 5, and in other embodiments, the multiplication factor a is substantially equal 
to 3. 

[0053] The acceleration sensors provide the following measurements at epoch k: 

A uk ] T ; (Eq. 13) 

where A Lf k is the measurement from the lower acceleration sensor 50 and Ay,k is the 
measurement from the upper acceleration sensor 40 at epoch k. The co-variance matrix 
corresponding to the acceleration sensor measurements at epoch k is provided by the 
following diagonal matrix: 



0 



(Eq. 14) 



where a A k is the uncertainty of the lower acceleration sensor measurements and a ^ k is the 
uncertainty of the upper acceleration sensor measurements. In certain embodiments, a A k is 
the same for all epochs, and is the same for all epochs. In certain embodiments in 
which the two acceleration sensors are substantially identical, cr A k = £7 a k . 
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[0054] As discussed more fully below, in certain embodiments, additional aiding 
data may be supplied. In certain embodiments, aiding velocity measurements [Vk] are 
provided with a corresponding co-variance matrix ^ fA =[o"£ t4 ], where <j Vk is the 
uncertainty of the aiding velocity measurements. In certain embodiments, aiding absolute- 
depth measurements [S*] are provided with a corresponding co-variance matrix i i J s k = [cr^ J, 

where cr S k is the uncertainty of the aiding absolute-depth measurements. In certain 
embodiments, aiding relative-depth measurements [Rk] are provided with a corresponding co- 
variance matrix 4^ =[c"^*]> where a R k is the uncertainty of the aiding relative-depth 

measurements. Other aiding measurements can be used in accordance with embodiments 
described herein. 

[0055] The theoretical acceleration sensor measurements can be calculated using 
the predicted state vector elements a* and h in the following equations: 

A \k = <*k + g * cos(7*) ; (Eq. 15) 

A'ujt= fl* + g*eos(/*-rf**£); (Eq. 16) 

where A ' Lik is the theoretical lower acceleration sensor measurement, A 'u,k is the theoretical 
upper acceleration sensor measurement. 

[0056] The equations which provide the predicted state vector at epoch k based on 
the state vector at epoch k-1 can be expressed as the following prediction matrix 



(Eq. 17) 



In this way, the predicted state vector at epoch k can be expressed as X k = Ok * Xk-h 

[0057] The prediction matrix & k can not be used for a prediction of the state co- 
variance matrix because it is non-linear (i.e., one of the state elements is included in the 
matrix). Instead, a linear prediction matrix 71 can be used to update the state co-variance 
matrix, corresponding to the uncertainty of the new state vector, as follows: 

I,=r i *E i _ 1 *r i r +^; (Eq.18) 
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where 7* is given by the following: 

1 0 0 0 0" 

At 10 0 0 

T k = {Atf/2 At 1 0 0 . (Eq. 19) 

0 0 0 1 0 

0 d k _i*At 0 v k _ } *At 1 

[0058] Three matrices can be defined to be used to calculate updates of the state 
vector Xk and the state co-variance matrix 2* based on measurements. The design matrix a A .k 
corresponds to the partial derivatives of the theoretical measurements and is given by the 
following: 

1 0 0 0 -g*sin(/ 4 _) 

10 0 sm{l k _-d k _*B)*B - sin(/ t . - d k _ * B) 

The constant vector fi A ,k corresponds to the measurements minus the theoretical 
measurements and is given by the following: 

- - ^*cos(/ jt _) 



a 



(Eq. 20) 



' Kk 






Kk - 


_A{j,k 






_ A u,k -0*- 



(Eq. 21) 



The gain matrix Gk is given by the following: 

G k = 2,. * a\ k * (v A<k + a A<k * I t _ * a \ k )"' . (Eq. 22) 

The after-measurement update of the state vector X k and the after-measurement update of the 
state co- variance matrix Z k are calculated as follows: 

X k = X k _ - G k * (a Atk * X k _ + p AJk ) ; (Eq. 23) 

JL k =* k --G k *a AJt *2 k _. (Eq.24) 

where k- denotes the values of epoch k prior to the current update. Using such a labeling 
scheme, k- in Equations 20 through 24 denotes the values of the predicted state vector prior 
to the update using the acceleration measurements. Thus, X k . and E k . are the state vector and 
state co-variance matrix, respectively, of epoch k after the prediction but prior to the 
measurement update. 
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Aiding Data 

Stationary Data 

[0059] In certain embodiments, the downhole portion 30 of the survey tool 10 is 
stopped within the wellbore 20 one or more times to obtain aiding data while the downhole 
portion 30 is stationary. In certain embodiments, the downhole portion 30 is stopped at one 
or more random or arbitrary times, while in other embodiments, the downhole portion 30 is 
stopped at multiple times with a predetermined period. In certain such embodiments which 
only have the first acceleration sensor 40 and the second acceleration sensor 50, the state 
vector Xk and the co-variance matrix can be expressed as: 

(Eq. 25) 
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which in certain embodiments is given by: 



(Eq. 26) 
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(Eq. 27) 



where k- in Equations 25 through 27 denotes the values of the various elements prior to the 
update using the aiding stationary data. In other such embodiments which have cross-axial 
acceleration sensors, the state vector Xk and the co-variance matrix 2* can be expressed using 
the following: 

d k = (arctan(- HA L k /A LJc )- arctan(- HA UJc fA UM ))/B ; (Eq. 28) 

I k =nctzn(-HA L jA Lk ); (Eq. 29) 
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2* = 



0 0 0 
0 0 0 

0 0 0 



0 
0 
0 



0 0 0 <r 2 ,j{Bji) 
which in certain embodiments is given by: 



0 
0 
0 



(Eq. 30) 
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(Eq.31) 



Velocity Data 

[0060] In certain embodiments, aiding velocity data is provided by one or more 
supplementary sensors which are part of the downhole portion 30 of the survey tool 10. As 
described above in relation to Figure 3, exemplary supplementary sensors include, but are not 
limited to, gamma-ray sensors adapted to detect gamma rays from geological formations in 
proximity to the downhole portion 30 and magnetic sensors adapted to detect casing collars 
in proximity to the downhole portion 30. Such supplementary sensors do not continuously 
measure the velocity of the survey tool 10, but they detect landmark features which can be 
used to calculate the average velocity of the downhole portion 30 as the downhole portion 30 
passes the landmark feature or between the locations of two or more selected landmark 
features. 

[0061] Figure 5 schematically illustrates the signals from a first gamma-ray sensor 
70 and a second gamma-ray sensor 80 as the downhole portion 30 moves past a geological 
formation 90 having a gamma-ray emission higher than the surrounding formations. The 
first gamma-ray sensor 80 is lower than the second gamma-ray sensor 70, so the first gamma- 
ray sensor 80 detects the formation 90 before the second gamma-ray sensor 70 detects the 
same formation 90. By correlating the measurements of the two supplementary sensors, and 
knowing the distance B v between the two supplementary sensors, the average velocity of the 
downhole portion 30 as it passed the geological formation 90 can be calculated as follows: 

V^BJ^-T,)- (Eq.32) 



-19- 



where B v is the distance between the first gamma-ray sensor 80 and the second gamma-ray 
sensor 70, T L is the time that the first gamma-ray sensor 80 detects the geological formation 
90, and Tu is the time that the second gamma-ray sensor 70 detects the geological formation 
90. Similarly, two magnetic sensors of the downhole portion 30 can be used to determine an 
average velocity of the downhole portion 30 relative to a casing collar of the wellbore 20. 

[0062] The measured velocity Vk is the average velocity over the time interval 
ranging from Ti to Tu. However, the average velocity approaches the instantaneous velocity 
as the time interval (T L -Tu) approaches zero. In certain embodiments, the average velocity is 
used to analyze the state vector X^ In certain other embodiments, the analysis of the state 
vector Xk utilizes the instantaneous velocity at epoch k. In certain such embodiments, the 
distance B v is selected to be as small as possible without letting the noise or uncertainty in the 
two time measurements corrupt the measured velocity. 

[0063] The uncertainty oy,* in the measured velocity at epoch k depends on the 
instability of the velocity over the time interval between T L and T u , and on the uncertainty in 
the repeated detection of the landmark feature being referenced. An estimate of the 
uncertainty of the measured velocity is given by: 



where Gdet.k is the uncertainty in the detection of the actual location of the formation at epoch 
k, and p v is the maximum change of the velocity of the time period At of the epoch k. In 
certain embodiments, <Jdet t k is constant among different epochs. In certain embodiments, the 
maximum change of the velocity over the time period At is assumed to be equal to the 
maximum change of the depth po divided by the time period At, where At is small. 

[0064] In certain embodiment in which such aiding velocity data is available, a 
design velocity vector ay t k (corresponding to the partial derivatives of the theoretical velocity 
measurements which are equal to the current state velocity) and a constant velocity vector 
pv.k (corresponding to the measurements minus the theoretical measurements) can be 
expressed as: 




(Eq. 33) 



a Vk =[0 1 0 0 0]; 



(Eq. 34) 



(Eq. 35) 
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Using these two vectors, a gain velocity matrix G* can be expressed as: 

G k = * a T rjk * (%,* + <v* * * <* I" 1 ; (Eq- 36 ) 

which can be used to express the state vector Xk and the co-variance matrix 2* as: 
X k = X k _ - G k * (a,,, * ^ + P v , ) ; (Eq. 37) 

Z 4 = Z,_-G,*a^*S A _; (Eq. 38) 

where £- in Equations 35 through 38 denotes the values of the various elements prior to the 
update using the aiding velocity data. 
Absolute Depth 

[0065] In certain embodiments in which the absolute depth of a landmark feature 
is known, the detection of this landmark feature by one or more sensors of the downhole 
portion 30 can be used to provide aiding data. In such embodiments, the time at which the 
landmark feature is passed by the downhole portion 30 is noted, and for the corresponding 
epoch the depth S* can be expressed as: 

S k =S s k +B s ; (Eq.39) 

where S% is the depth of the landmark feature and B s is the distance between the lower 

acceleration sensor 50 and the supplementary sensor in proximity to the lower acceleration 
sensor 50 (e.g., the second gamma-ray sensor 80 as in Figure 3). 

[0066] In certain embodiments, one or more of the same supplementary sensors 
are used to provide both aiding velocity data and aiding absolute-depth data. In certain 
embodiments, two supplementary sensors are used to provide two separate absolute-depth 
measurements of the same landmark feature. In such embodiments, the two measured 
absolute depths, S* and are related to two different epochs which are temporally close, 
and are related by the following: 

S k = S s k +B SL ; (Eq.40) 

S k ,=S s k +B su ; (Eq.41) 

where B SL is the distance between the lower acceleration sensor 50 and the second 
supplementary sensor (e.g., the second gamma-ray sensor 80 as in Figure 3), and B? u is the 
distance between the lower acceleration sensor 50 and the first supplementary sensor (e.g., 
the first gamma-ray sensor 70 as in Figure 3). 
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[0067] The uncertainty as , k of the absolute-depth measurements at epochs k and k 9 
depends on both the uncertainty in the given value of the absolute depth of the landmark 
feature and the uncertainty in the detection of the landmark feature by the supplementary 
sensor. The uncertainty of the absolute-depth measurement can be expressed as follows: 

<* =^dea + ^v^; (Eq.42) 
where a detlk is the uncertainty in the detection of the landmark feature and a^enM is the 
uncertainty in the given value of the absolute depth of the landmark feature. In certain 
embodiments, a^k is constant for different epochs. 

[0068] In certain embodiments in which aiding absolute-depth data is available, a 
design absolute-depth vector a s ,k, a constant absolute-depth vector fis.h and a gain absolute- 
depth matrix G k can be used to express the state vector X k and the co-variance matrix E k as 



follows: 

a S k =[0 0 1 0 0]; (Eq. 43) 

Ps*AS k -D k _\ (Eq.44) 

G k = E A _ * <* * K* + a SJc * E t _ * a T SJk Y ; (Eq. 45) 

X k =X k _-G k * (a SJt * X k _ + fi 8Jt ) ; (Eq. 46) 

Z 4 =I,.-G,*a^*2,_; (Eq.47) 



where in Equations 44 through 47 denotes the values of the various elements prior to the 
update using the aiding absolute-depth data. 
Relative Depth 

[0069] In certain embodiments, the downhole portion 30 provides relative-depth 
data which does not depend on the known depths of landmark features. In certain 
embodiments, such relative-depth data is provided by two supplementary sensors adapted to 
provide large and sharp signal spikes at landmark locations. In certain embodiments, these 
landmark locations are previously known, while in other embodiments, these landmark 
locations are previously unknown. In certain embodiments, the relative-depth data is 
provided by a similar configuration of two supplementary sensors of the downhole portion 30 
as that which provides the aiding velocity data described above. However, while the velocity 
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data is advantageously provided by two supplementary sensors which are relatively close 
together, the relative-depth data is advantageously provided by two supplementary sensors 
which are farther apart from one another. 

[0070] In certain embodiments, the epoch k is defined to be the epoch during 
which the second supplementary sensor (e.g., the second gamma-ray sensor 80 of Figure 3) 
passes a detectable landmark location and the epoch k } is defined to be the epoch during 
which the first supplementary sensor (e.g., the first gamma-ray sensor 70 of Figure 3) passes 
the same detectable landmark location. The aiding relative-depth measurement Rk at epoch k 
of certain embodiments can be expressed as: 



where Z)*> is the calculated depth of the downhole portion 30 at epoch k\ B r is the distance 
between the first supplementary sensor and the second supplementary sensor (e.g., the two 
gamma-ray sensors 70, 80 of Figure 3), and B a is the distance between the second 
acceleration sensor 50 and the second supplementary sensor (e.g., the second gamma-ray 
sensor 80 of Figure 3). 

[0071] In certain embodiments, the uncertainty in the relative-depth measurement 
at epoch k depends on both the uncertainty of the calculated depth at epoch k' and the 
uncertainty in the relative depth detection. The uncertainty in the relative-depth measurement 
can be expressed as follows: 



where <j re i is the uncertainty in the relative-depth detection for the landmark location and 
(7iast,k is the uncertainty of the calculated depth at epoch k. In certain embodiments, <Ji as a is 
given by the third diagonal element (2jj f f) of the solution co-variance matrix at epoch In 
certain embodiments, oiasu is equal to <jd, which is the uncertainty in the depth measurement. 

[0072] In certain embodiments in which aiding relative depth data is available, a 
design relative-depth vector ccs.h a constant relative-depth vector and a gain relative- 
depth matrix G* can be used to express the state vector Xk and the co-variance matrix 2J as 
follows: 



R k =D k ,+B r +5 fl ; 



(Eq. 48) 




last,k > 



(Eq. 49) 



a 



w=[0 0 1 0 0]; 



(Eq. 50) 
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t 



(Eq.51) 



(Eq. 52) 



(Eq. 53) 



2* =z*--G 4 s*.; 



(Eq. 54) 



where £- in Equations 51 through 54 denotes the values of the various elements prior to the 
update using the aiding relative-depth data. 

[0073] Various embodiments of the present invention have been described above. 
Although this invention has been described with reference to these specific embodiments, the 
descriptions are intended to be illustrative of the invention and are not intended to be limiting. 
Various modifications and applications may occur to those skilled in the art without departing 
from the true spirit and scope of the invention as defined in the appended claims. 
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